ID-A133  089  THE  EFFECT  OF  FIELD  INDEX  FLUCTUATIONS  ON  THE  DVNflMICS 

OF  A  HIGH  CURRENT.  .  <U)  NAVAL  RESEARCH  LAB  WASHINGTON  DC 
C  AGRITELLIS  ET  AL.  14  SEP  82  NRL-MR-5141 


UNCLASSIFIED 


F/G  9/2 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


4  TjTuE  'end  Subtitle) 

THE  EFFECT  OF  FIELD  INDEX  FLUCTUATIONS  ON  THE 
DYNAMICS  OF  A  HIGH  CURRENT  ELECTRON  RING  IN 

A  MODIFIED  BETATRON  ACCELERATOR 

i  type  of  report  <  PERIOD  COVERED 

Interim  report  on  a  continuing 

NRL  problem. 

«  PERFORMING  org-  report  number 

7  AUTHOR' »/ 

C.  Agriteilis,*  S.J.  Marsh, ♦♦  and  C.A.  Kapetanakos 

8.  contract  or  grant  number^; 

t  PERFORMING  ORSAni  Z  *YiOn  NAME  INO  ADDRESS 

Naval  Research  Laboratory 

Washington,  DC  20375 

•o  program  element,  project,  task 
AREA  &  WORK  UNIT  NUMBERS 

61153N;  RR011-09-41; 

47-1484-0-3 

1  1.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Office  of  Naval  Research 

12.  REPORT  0 ate 

September  14,  1983 

Arlington,  VA  22217 

13.  NUMBER  OF  PAGES 

30 

l«.  monitoring  AGENCY  name  4  AOORESS (II  Otlloront  Iron  Controlling  Olllco) 

19.  5CC  U  R  t  T  Y  CLASS.  t ot  thle  report) 

UNCLASSIFIED 

1S«.  DECLASSIFICATION'  DOWNGRADING 

schedule 

_ 1 

i«.  distribution  statement  rot  tnt a  R*ponj 

Approved  for  public  release;  distribution  unlimited. 


>7.  DISTRIBUTION  STATEMENT  ;ot  the  abatrect  entered  In  Block  20,  It  different  (toot  Report) 


it.  supplementary  notes 

♦Present  address:  Science  Applications,  Inc.,  McLean,  VA  22102 
♦♦Present  address;  Sachs/Freeman  Associates,  Bowie,  MD  20715 


19  KEY  WOROS  (Continue  on  rover  am  •  ‘dm  ti  necee  aery  mnd  Identify  lay  block  number) 

Modified  betatron 

Betatron 

Accelerators 

High  current  accelerators 


20.  ABSTRACT  ( Continue  on  reveree  aide  it  neceeeery  and  identity  by  block  number) 

y  Analytical  and  computer  simulation  results  are  reported  on  the  effect  of  the  external  field  index 
fluctuations  on  the  dynamics  of  a  high  current  electron  ring  in  a  modified  betatron.  These  non-linear 
studies  clearly  demonstrate  that  the  external  field  index  fluctuations  are,  in  general,  harmless  to  the 
high  current  ring.  In  addition,  it  was  found  that  a  small  (1%)  temperature  spread  parallel  to  the 
direction  of  propagation  excites  oscillations  in  the  rros  emittance,  but  the  amplitude  of  these 
oscillations  is  substantially  lower  than  in  a  conventional  betatron. 


FORM 

I  JAW  71 


1473  COlTlON  OF  I  NOV  «9  IS  obsolete 

S/N  0  102-0  14-  <40  1 


SECURITY  CLASSIFICATION  OF  THIS  RASE  (Whtn  Omlm  tn  lotoel 


THE  EFFECT  OF  FIELD  INDEX  FLUCTUATIONS  ON  THE 
DYNAMICS  OF  A  HIGH  CURRENT  ELECTRON  RING  IN  A 
MODIFIED  BETATRON  ACCELERATOR 


I.  Introduction 

The  linear  dynamics  of  a  high  current  electron  ring  in  a  modified  betatron 

geometry  has  been  studied  extensively*”**  during  the  last  few  years.  These 

studies  are  mainly  motivated  by  extensive  evidence  suggesting  that  the 

modified  betatron  has  the  potential  to  confine  high  current  electron  rings, 

12 

*eyen  in  the  presence  of  a  substantial  transverse  emittance 

The  modified  betatron  belongs  to  the  class  of  cyclic  induction 
accelerators,  i.e.,  the  electrons  are  accelerating  by  an  inductive  electric 
field.  So  far  these  different  cyclic  induction  accelerators  have  been 
proposed:  the  conventional  betatron*3”*5,  the  modified  betatron*”**  and  the 

stellatron*5.  The  modified  betatron  includes  in  addition  to  the  time  varying 
betatron  magnetic  field  that  is  responsible  for  the  acceleration,  a  strong 
toroidal  magnetic  field  that  substantially  improves  the  stability  of  the 
conventional  betatron.  In  the  stellatron,  the  addition  of  a  stellarator  field 
to  the  modified  betatron  substantially  reduces  the  displacement  of  the  orbit 
that  is  due  to  energy  mismatch.  However,  beam  trapping  and  resonances  appear 
to  be  major  and  presently  unsolved  problems. 

The  non-linear  dynamics  of  very  tenuous  electron  rings  in  a  modified 
betatron  accelerator  has  been  studied  recently  by  Chernin*^.  His  studies  were 
limited  to  quadratic  non-linearities  in  the  absence  of  self-fields  and  image 
fields  from  the  surrounding  walls.  For  azimuthally  symmetric  fields,  he 
concluded  that  the  quadratic  non-linearities  are  quite  harmless  to  the 
electron  ring.  In  this  paper,  we  present  computer  simulation  results  which 

clearly  demostrate  that  Chernin's  conclusion  is  valid  under  much  more  general 
Manuacript  approved  June  28,  1983. 
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conditions,  i.e.,  for  high  current  rings  in  the  present  of  surrounding  walls 
and  strong  non-linearities  as  for  example  when  the  spatial  derivatives  of  the 
external  field  index  are  large. 


II.  Computer  Simulation  Results 


a.  Description  of  the  particle  in  cell  computer  code. 


The  computer  code  has  been  described  previously0.  Here  we  provide  a  very 
brief  summary  with  emphasis  on  a  recent  modification  related  to  the  initial 
loading  of  the  particles.  Briefly,  the  particle  simulation  code  is  2D  in 
configuration  and  3D  in  velocity  space.  It  computes  self-consistently  all 
self  field,  except  the  self  Bg  field.  This  assumption  is  valid  provided 
v/y  is  small.  In  addition  the  radiative  term  (displacement  current)  is 
ignored,  i.e.,  the  code  uses  the  Darwin  model  for  Maxwell’s  equations. 

The  electrostatic  potential  is  computed  from  Poission’s  equation 


v  ♦  ■  p/e  , 
.  o 


and  the  magnetic  vector  potential  from 
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with  the  boundary  condition  $  *  Aq  ■  0  at  the  conducting  wall. 
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Equations  (1)  and  (2)  are  solved  by  Fourier  decomposition  in  the  z- 

direction  and  then  by  Gaussian  elimination  of  the  resultant  tridiagonal  matrix 

2 

of  equations  obtained  from  a  3-point  differencing  scheme  for  7^  .  The  inverse 
Fourier  transform  yields  Ag  and  $  on  the  grid.  The  particle  velocities  in  the 
toroidal  direction  are  obtained  using  the  conservation  of  cannonical 
momentum.  Therefore,  the  equation  for  Ag  is  not  properly  time  centered  since 
the  velocities  from  the  previous  time  step  are  used  to  calculate  the  currents 
from  the  cannonical  momenta.  This  method  was  chosen  primarily  for  its  speed 
and  simplicity  but  care  must  be  taken  in  applying  the  code  when  the  inductive 
acceleration  of  particles  in  the  toroidal  direction  is  significant. 

In  several  runs  the  particles  were  loaded  in  the  code  at  t  *  o  using  a 
cylindrical  K-V  distribution.  In  these  runs  the  electrons  quickly  acquired  an 
“energy  spread".  This  “thermalization"  is  due  to  the  fact  a  cylindrical  K-V 
distribution  is  not  appropriate  for  high  current  electron  rings  that  have 
large  aspect  ratio  rQ/rb,  where  rQ  is  the  major  and  rb  the  minor  radii  of  the 
ring.  The  reason  is  that  a  cylindrical  K-V  distribution  requires  that  the 
potential  be  the  same  on  both  the  left  and  right  side  of  the  beam.  However, 
in  a  actual  ring  there  is  a  potential  difference  between  the  inner  and  outer 
edge.  Thus,  a  ring  that  has  been  incorrectly  initialized  tries  to  attain  a 
more  physical  distribution,  but  in  the  absence  of  dissipation  this  can  be 
achieved  only  temporarily.  In  the  process  a  spread  in  y  is  developed,  which 
is  equivalent  to  temperature.  In  all  the  results  presented  in  this  paper  the 
code  was  loaded  at  t  -  o  correctly,  i.e.,  taking  into  account  the  asymmetry  of 
the  potential  across  the  ring. 


b.  Results  without  temperature  in  the  azimuthal  direction. 

Figure  1  shows  the  modified  betatron  geometry  and  the  system  of 

coordinates  that  will  be  used  in  our  subsequent  discussion.  The  center  of  the 

electron  ring  is  located  away  from  the  center  of  the  minor  cross-section  of 
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the  torus  and  at  a  distance  A  *  (Ar  +  Az  )  .  Results  from  the  computer 

simulation  are  shown  in  Fig.  2.  The  various  parameters  for  this  run  are 
listed  in  Table  I.  At  t  *  o  the  5  KA  electron  ring  is  located  on  the 
horizontal  symmetry  plane  and  at  a  distance  A  *  8  cm,  i.e.,  half-way  between 
the  minor  axis  and  the  wall  of  the  toroidal  chamber.  It  is  apparent  from  the 
small  variations  of  the  rings  envelopes  l Figs.  2a  and  2b]  that  the  ring  is 
reasonably  well  matched.  As  shown  in  Table  I  the  external  field  index  n  is 
constant  and  equal  to  0.42.  For  such  a  value  of  n  the  linear  theory  predicts 
that  the  minor  axis  of  the  ring  has  two  modes  of  oscillation,  a  slow  (bounce 
frequency)  and  a  fast  (cyclotron  frequency  corresponding  to  the  toroidal 
field).  The  orbit  of  the  ring  axis  associated  with  the  slow  or  bounce  mode  is 
approximately  a  circle  that  is  centered  around  the  ring's  equilibrium 
position.  Although  the  run  was  terminated  before  the  ring  could  complete  a 
full  bounce  period,  it  is  clear  from  Fig.  2c  that  the  center  of  the  ring 
describes  an  orbit  that  is  similar  to  that  predicted  by  the  theory.  Both 
modes  of  oscillation  are  quite  apparent.  The  radius  of  the  fast  mode  can 
be  computed  approximately  from  the  expression  p^  *  (v^c^/r^Jlg) ,  w^ere  ve  *s 
the  azimuthal  velocity,  is  the  normalized  emittance,  r^  is  the  ring  radius 


and  S2g  the  cyclotron  frequency  of  the  local  toroidal  magnetic  field.  The  time 
interval  between  two  successive  arrows  in  the  figure  is  20  nsec. 

The  ring  kinetic  energy  also  oscillates  with  the  two  characteristic  modes 
as  may  be  seen  in  Fig.  2d.  The  variations  of  y  can  be  computed  approximately 
by  integrating  the  energy  rate  equation 


dy  *  *  0 

3t  ’  '  el  V-E  ’  el  VH’ 


using  for  the  self  magnetic  vector  potential  the  expression 


A0  *  ~c~  t1/2  +  ln  (a/rb)  +  ltt  *  A2/ a2) ) • 


In  Eq.  (4),  I  is  the  ring  current,  a  is  the  minor  radius  of  the  torus 
and  A  is  the  distance  between  the  ring's  minor  axis  and  that  of  the  torus.  It 
should  be  noticed  that  Eq.  (4)  does  not  include  toroidal  effects  and  therefore 
is  valid  only  for  low  v/y. 


Since  If  ■  -  |e|  N  Vg^nr'  w^ere  ^  tota*  number  of  electrons  in  the 

ring  and  the  external  part  of  the  magnetic  vector  potential  remains  constant, 
we  get  from  Eqs.  (3)  and  (4) 


y  (t)  -  y  (o)  *  2vq  {-p™-  (1/2  +  in  (a/rfe  (o))+ 


in  (1  -  A2(o)/a2)]  -  -A-  (1/2  +  in  (a/r.  (t))  +  in  (1  -  A2(t)/_2)  ]} ,  (5) 


where  vq  is  the  Budker  parameter  at  rQ,  i.e., 


In  the  fast  mode  y  changes  through  the  transverse  oscillations  of  the  ring's 
envelope.  Since  in  this  fast  time  scale  r  and  A  remain  approximately 
constant,  Eq.  (5)  becomes 

Y  (t)  -  y  (o)  -  2v  —■  in  [r.(t)/r  (o)J. 

or  D  D 

/ 

In  the  slow  mode  y  changes  mainly  through  the  variation  of  the  ring 
current,  which  results  from  the  changing  r.  Since  the  contribution  from  the 
variation  of  r^  and  A  are  in  general  small,  Eq.  (5)  gives 

y  (t)  -  y  (o)  «  2vq  rQ  [1/2  +  in  (a/rb(o>)  + 

ta  (1  -42(.)/,2)l  (7) 

Equation  (7)  predicts  that  y  (t)  decreases  with  r  (t).  In  addition,  for 
the  results  of  Fig.  2d,  Eq.  (7)  predicts  a  reduction  in  y  between  t  *  o  and 
t  »  80  nsec  of  0.18,  which  is  in  excellent  agreement  with  the  computer 
simulation  results. 

As  it  was  stated  before,  the  run  shown  in  Fig.  2  was  made  using  a  constant 
external  field  index.  The  purpose  this  run  was  to  serve  as  a  "bench  mark" 
for  the  rest  of  the  runs  that  were  made  with  a  variable  external  field 


index.  The  radial  profile  of  the  betatron  field  and  the  variation  of  the 
external  field  index  are  shown  in  Fig.  3.  The  assumed  variation  is  in  effect 
greater  than  that  expected  to  be  present  in  a  well  designed  device. 

The  effect  of  the  external  field  index  variation  on  the  dynamics  of  the 
ring  is  shown  in  Fig.  A.  The  parameters  of  this  run  are  listed  in  Table  II 
and  with  the  exception  of  the  external  field  index  are  identical  to  those  of 
run  listed  in  Table  I.  By  comparing  Figs.  (2)  and  (A)  it  may  be  concluded 
that  the  variation  of  the  field  index  does  not  have  a  profound  effect  on  the 
dynamics  of  the  ring.  Although  some  details  are  different  the  gross  features 
of  the  two  runs  are  very  similar.  The  most  pronounced  new  feature  of  the 
results  in  Fig.  A  is  the  slow  time  scale  variation  of  the  ring  envelope.  As 
the  ring  moves  from  its  initial  position  at  z  ■  o  to  z  *  o,  the  field  index  is 
reduced  and  the  ring  becomes  unmatched,  resulting  in  envelope  oscillations. 

Snap-shots  of  the  minor  cross  section  of  the  ring  and  the  magnetic  field 
lines  corresponding  to  the  total  magnetic  field  are  shown  in  Fig.  5.  As  can 
be  easily  computed,  the  zero  magnetic  field  point  occurs  at  ~  rb/2 

Similar  results  to  chose  shown  in  Fig.  A  were  also  obtained  for  smaller 
initial  ring  displacements  A  (o).  Figure  6  shows  the  results 
for  A  (o)  ■  5  cm.  The  various  parameter  of  this  run  are  listed  in  Table 
III.  In  all  the  runs  we  carefully  avoided  crossing  the  beam  resonance,  i.e., 
the  radial  frequency 

V  (T£)2  t^2-n  (t)  5-  Vb2/a2], 


where  i  ■  U  +  2  (v/^,)  (0.5  +  in  a/r^)] 

was  kept  different  than  zero.  The  implications  of  crossing  the  resonance  are 
presently  under  investigation. 

-o.  The  effect  of  temperature 

In  the  system  of  coordinates  shown  in  Fig.  1,  the  equations  describing  the 

motion  of  individual  electrons  in  a  constant  radius  beam  with  A  ■*  o  having  an 

azimuthal  energy  spread  Ay/  are 

Yo 

2 

6r  +  J  5r  -  <0,9/ V  *  *  *  7-  (8) 

o  To 

5z  +  w*  <Sz  +  <flo0/Yo)  Sr  «  o,  (9) 

2  2  2  2 
where  w  «  (fl  /y  )  (1-n-n  ) ,  u  -  (Q  /y  )  (n-n  )  and  Ay  =  Y~  <y>  is  the 

r  oz  o  s  z  oz  o  s 

azimuthal  energy  spread  in  the  beam.  Equations  (8)  and  (9)  do  not  include  the 
toroidal  corrections  in  the  fields  and  therefore  are  valid  only  for  low 
v/Y0  beams. 

In  the  Larmor  frame,  i.e.,  a  frame  that  rotates  with  a  constant  frequency 

Q,  *  1/2  (Q  J Y  )»  Eqs.  (8)  and  (9)  take  the  form 
L  00  0  2  Rs 

«r+  .  J-  41  co.  (->4. 

r0V8Z  °  9 

A  y  >  A  A  2  .  S 

6z  +  K  6z  *  - =•  —  sin  ( - ), 


(11) 


<$z  ■  sin  ( - )  <Sr  +  cos  ( - )  $z. 


K2  -  ( oj2  +  Q2)  /  vj!,  ui2  *  a)  2  »  to  2,  i.e.,  n-1/2,  6r"  -  s  -  r  0  and  9  is 

L  o  r  z  as  o 

the  toroidal  angle. 

Equations  (10)  and  (11)  describe  forced  oscillations  and  have  a  resonance 
at  (D  =*  o. 

The  solution  of  Eqs.  (10)  and  (11),  when  K  is  independent  of  s  is 


6r  *  [6r(o)  -  5r  J  cos  Ks  +  sin  Ks  +  Sr  cos  (— — ) ,  (12) 

OK  O  V . 


*  *  1  A  ^  Qj  S 

6z  »  <Sz  (o)  cos  Ks  +  rr  [  6z  (o)  -  <5r  ft_/vQ]  sin  Ks  +  5r  sin  ( - ),  (13) 

K  o  L  0  o  Vq 


where 


2  . 
c  at 

:o  *  2  Y  * 

r  u>  'o 
o 


The  rms  beam  emittance  is  defined  by 

2 

e  -  4  {<5r  (t)>  <5r  (t)>  -  <6r(t)  5r  (t)»A/*. 


(14) 


Substituting  Eqs.  (12)  and  (13)  into  Eq.  (14)  and  taking  averages  over  a  K-V 
distribution,  we  obtain 


—  =■  <dr2(o)>  <ir'  (o)> 


■+5r  2{<5r2(o)>[ xKsinKs+x  cosKs]2-K6r  (o)>[xcosKs-£  sin  Ks]2}, 

O  Iv 


where 


ft^s 

cos  ( - )  -  cos  Ks 

v0 


nL  nT  s 

—  sin  ( - )  +  KsinKs, 

v0  v9 


for  the  r-component 


and 


V 

sin  (^-) 
0 


,QL.  $JLs 

( — )  cos  - 

v9  V0 


(-jT— ■)  sinKs 
Kv0 


(— )  cosKs 
V0 


for  the  z-component 


If  the  initial  conditions  are  such  that 


K2  <<5r2(o)>  -  <6r"2(o)>. 


(15) 


(16) 


(17) 


Eq.  (15)  becomes 


(18) 


Results  from  the  computer  simulation  with  a  small  parallel  temperature 
spread  are  shown  in  Fig.  7.  The  various  parameters  for  this  run  are  listed  in 
Table  IV.  When  a  parallel  energy  spread  is  present,  the  matching  becomes  more 
difficult  and  the  beam  envelope  starts  to  oscillate.  These  oscillations 
result  in  the  growth  of  the  ring's  rms  emittance  and  is  discussed  later  on. 
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The  ring  dynamics  in  the  presence  of  temperature  is  similar  to  that  in  the 

absence  of  temperature,  provided  that  even  with  the  expansion  of  the  beam 
2 

radius  oj  remains  different  than  zero.  In  one  run  the  various  parameters  were 

2 

assigned  such  values  that  w  would  go  through  zero  after  a  small  expansion  of 
the  beam  radius.  A  very  rapid  increase  in  the  beam  envelope  was  observed  that 
-was  accompanied  by  a  15-fold  growth  of  the  rms  emittance  in  approximately  100 
nsec.  However,  the  expansion  had  no  noticable  effect  on  the  macroscopic 
motion  of  the  center  of  the  beam,  at  least  for  the  duration  of  the  run. 

The  rms  emittance  as  a  function  of  time  for  T^  t  o  is  shown  in  Fig.  8a. 
The  various  parameters  are  listed  in  Table  IV.  The  corresponding  result 
for  T^  *  o  is  shown  in  Fig.  8b.  With  the  exception  of  the  temperature,  the 
various  parameters  for  this  run  are  identical  to  those  listed  in  Table  IV. 
Since  in  the  numerical  results 

<5r2(o)>  -  <Sz2(o», 

and 

<dr2(o)>  -  <Sz"2(o)>, 

Equations  (19)  and  (20)  predict  that  when  »  tu 


i.e.,  e  and  e  are  either  in  phase  or  approximately  180°  out  of  phase.  Fig. 
r  z 

8a  shows  that  e  and  £  are  180°  out  of  phase,  although  the  condition 
r  z 
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&og  >>  u  is  not  satisfied.  This  appears  to  be  a  general  results  and  has  been 
observed  in  all  the  runs  with  temperature  or  equivalent  temperature. 

The  secular  emittance  growth  in  Fig.  8a  is  attributed  to  the  radial 
expansion  of  the  ring.  This  is  also  a  general  result  and  is  observed  in  all 
unmatched  rings  with  or  without  temperature.  Among  the  various  non- 
linearities,  such  as  field  fluctuations,  wall,  density  and  radial  expansion, 
the  latter  appears  to  be  the  most  pronounced. 
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III.  Conclusions 


18 

Some  very  preliminary  experimental  results  together  with  the  extensive 
2-11  19-24 

equilibrium  and  stability  analysis  clearly  indicate  that  the  modified 

betatron  concept  has  the  potential  to  lead  to  the  development  of  very  high 
Current  accelerators. 

Our  confidence  in  the  modified  betatron  concept  is  further  enhanced  by  the 
present  results  which  show  that  the  magnetic  field  fluctuations,  that  are 
inevitable  in  any  practical  device,  appear  to  be  harmless  to  the  high  current 
ring.  However,  the  resonance  that  occurs  when  u  *  o  is  presently  a  major 
concern. 
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Table  I 


Parameters  of  the  run  shown  in  Fig.  2 

Run  No.  M14D 

Initial  beam  energy  yq  ■  7.117 

Beam  current  I  (KA)  *  5  KA 

Major  radius  rQ  (cm)  *  100 

Initial  beam  minor  radius  r^  (cm)  *  3 

Torus  minor  radius  a  (cm)  ■  16 

Initial  beam  center  position  r^  (cm)  *  108 

Betatron  magn.  field  at  rQ,  z  *  o,  Boz  (G)  ■  136.2 

Toroidal  magn.  field  at  rQ,  z  ■  o,  BoQ  (G)  »  388 

Initial  emittance  e  (rad  -  cm)  ■  0.1 

Initial  temperature  spread  (half-width)  ~  ■  o 

External  field  index  n  ■  0.42  (Constant) 

Self  field  index  n_  "  1.34 

9 

Wall  resistivity  p  *  0 
Time  step  (nsec)  *  100  ps 
No.  of  particles  -  1024 


Table  II 


Parameters  of  the  run  shown  in  Fig.  4 

Run  No.  Ml 6 A 

Initial  beam  energy  yq  *  7.117 

Beam  current  I  (KA)  -  5 

Major  radius  rQ  (cm)  *  100 

Initial  beam  minor  radius  r^  (cm)  *  3 

Torus  minor  radius  a  (cm)  ■  16 

Initial  beam  center  position  rA  (cm)  -  108 

Betatron  magn.  field  at  rD,  z  *  o,  BQ2  (G)  *  136.2 

Toroidal  magnetic  field  at  rQ,  z  -  o,  BqQ  (G)  *  388 

Initial  emittance  c  (rad  -  cm)  »  0.1 

Initial  temperature  spread  (half -width)  *  0 

External  field  index, n  *  see  Fig.  3 

Self  field  index  ng  *  1.34 

Wall  resistivity  p  ■  0 

Time  step  (nsec)  *  100 

No.  of  particles  *  1024 


Table  III 


Parameters  of  the  run  shown  in  Fig.  5 

Run  No.  M05A 

Initial  beam  energy  yQ  *  6.994 

Beam  current  I  (KA)  ■  4.76 

Major  radius  rQ  (cm)  *  100 

Initial  beam  minor  radius  r^  (cm)  »  3 

Torus  minor  radius  a  (cm)  *  16 

Initial  beam  center  position  r^  (cm)  -  105 

Betatron  magn.  field  at  rQ,  2*0  BQ2  (G)  *  136.2 

Toroidal  magn.  field  at  r0,  z  ■  o,  BqQ  (G)  *  388 

Initial  emittance  e  (rad  -  cm)  ■  0.1 

Initial  temperature  spread  (half-width)  ~  *  0 

External  field  index  n  ■  see  Fig.  3 

Self  field  index  ng  *  1.33 

Wall  resistivity  p  ■  0 

Time  step  (nsec)  *  100  ps 

No.  of  particles  ■  1024 
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Table  IV 


Parameters  of  the  run  shown  in  Fig.  7 

Run  No.  J02A 

Initial  beam  energy  yq  *  7.004 

Beam  current  I  (KA)  ■  5  KA 

Major  radius  rQ  (cm)  *  100 

Initial  beam  minor  radius  r^  (cm)  9  2 

Torus  minor  radius  a  (cm)  *  16 

Initial  beam  center  position  r^  (cm)  *  105 

Betatron  magn.  field  at  rQ,  z  »  o,  Boz  (G)  »  140 

Toroidal  magn.  field  at  rQ,  z  ■  o,  Bqq  (G)  -  725 

Initial  emlttance  e  (rad  -  cm)  -  0.1 

Initial  temerpature  spread  (half-width)  ~  ■  1Z 

External  field  index  n  ■  0.42 

Self  field  index  n_  ■  2.984 

Wall  resistivity  p  *  0 

Time  step  (nsec)  ■  100  ps 

No.  of  particles  *  1024 
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Fig.  2  (a)  Temporal  variation  of  rms  radial  distance  of  electrons  from  the 

center  of  the  minor  cross-section  of  the  ring;  (b)  temporal  variation 
of  the  rms  vertical  distance  of  electrons  from  the  center  of  the  minor 
cross-section  of  the  ring;  (c)  orbit  of  the  ring’s  center  in  the  r-z 
plane;  (d)  temporal  variation  of  ring's  kinetic  energy.  In  this  run 
the  external  field  index  is  uniform.  The  various  parameters  are 
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Variation  of  the  betatron  magnetic  field  and  external  field  index  with 


Root  tloon  8r 
S  I.SO  I.3t 


Fig.  5  Snap-shots  of  the  electron  ring  minor  cross-section  and  magnetic  field 
lines  corresponding  to  the  total  magnetic  field.  The  infinite 
conductivity  vacuum  chamber  has  a  circular  minor  cross-section.  The 
major  axis  of  the  torus  is  located  to  the  left  and  at  100  cm  from  the 
center  of  the  minor  cross-section  of  the  torus.  The  various 
parameters  are  listed  in  Table  II. 


Fig.  6  (a)  Temporal  variation  of  the  nns  radial  distance  of  electrons  from 

the  center  of  the  minor  cross-section  of  the  ring;  (b)  temporal 
variation  of  the  rms  vertical  distance  of  electrons  from  the  center  of 
the  minor  cross-section  of  the  ring;  (c)  orbit  of  the  ring's  center  in 
the  r-z  plane;  (d)  temporal  variation  of  ring's  kinetic  energy.  In 
this  run  the  external  field  index  varies  as  shown  in  Fig.  3.  The 
parameters  for  this  run  are  listed  in  Table  III. 


Z  leal 


Emittance  vs  Time 


Emittance  vs  Time 


Fig.  8  (a)  Temporal  variation  of  the  rms  emittance  e  ,  e  and  their  average 

r  z 

A  A 

(e  +  e  )/2,  when  there  is  a  1%  (half-width)  temperature  spread  T,. 
in  the  toroidal  direction.  The  various  parameters  for  this  run  are 
listed  in  Table  IV;  (b)  as  in  (a)  but  with  T, ,  ■  o. 
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